epi draft 



Enhanced backscattering of a dilute Bose-Einstein condensate 



GUILLAUME LABEYRIE^ , TOMASZ KARPIUK^'^ , 

TiAN MiNiATURA^'^'^ and Dominique Delande^ 



Jean-Francois Schaff^'^ , Benoit Gremaud^'^'^ , Chris- 



^ ^ Institut Non Lineaire de Nice, UMR 7335, UNS, CNRS; 1361 route des Lucioles, F- 06560 Valhonne, France 

Center for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, 
Singapore 



o 

(N 



a 

O 

o 



Wydzial Fizyki, Uniwersytet w Bialymstoku, ul. Lipowa 4^1, 15-424 Bialystok, Poland 
Vienna Center for Quantum Science and Technology, Atominstitut, TU Wien, 1020 Vienna, Austria 
Department of Physics, National University of Singapore, 2 Science Drive 3, Singapore 117543, Singapore 
Laboratoire Kastler Brossel, Ecole Normale Superieure, CNRS, UPMC; 4 Place Jussieu, F-75005 Paris, France 



PACS 67 . 85 . Hj - Bose-Einstein condensates in optical potentials 
PACS 42 . 25 . Dd - Wave propagation in random media 
PACS 05 . 60 . Gg - Quantum transport 

Abstract -We study experimentally and numerically the quasi-bidimensional transport of a *^Rb 
Bose-Einstein condensate launched with a velocity vo inside a disordered optical potential created 
by a speckle pattern. A time-of-flight analysis reveals a pronounced enhanced density peak in the 
backscattering direction —vo, a feature reminiscent of coherent backscattering. Detailed numerical 
simulations indicate however that other effects also contribute to this enhancement, including a 
"backscattering echo" due to the position-momentum correlations of the initial wave packet. 
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Introduction. — Wave localization in non- 
homogeneous media is a general phenomenon, based 
on interference of multiple partial waves. While it is well 
understood in spatially periodic media, with e.g. band 
gap properties for electronic waves in solids and photons 
in photonic crystals, its generalization to disordered 
media has sparked a wide interest since the pioneering 
work of Anderson [Ij in 1958. Anderson localization 
has been investigated in a variety of systems including 
light [2], microwaves [31 or ultrasound j^. Recently, 
dramatic progresses in the cooling and control of neutral 
atoms have allowed the investigation of the localization of 
quantum waves by disorder pjQl (see p!Ol for a review). 
Observing Anderson localization in dimension 2 or 3 
usually requires strong scattering, a condition difficult to 
realize experimentally. However, precursors of localization 
already exist in the regime of "weak localization" [11]. 
Coherent backscattering (CBS), a constructive interfer- 
ence effect in the backscattering direction, is a hallmark 
of this regime and can be observed in a wide range of 
transport mean free path i. It is a signature of (quan- 
tum) coherent transport and consequently an important 
tool for characterizing scattering media, especially for 



studies of localization phenomena. CBS has been studied 
extensively with light [12^ 13j, ultrasound [H], seismic 
waves [15], excitons [6] and polaritons [TT] . 

In this paper we report our experimental study of the 
transport of a Bose-Einstein condensate (BEC) in a bidi- 
mensional (2D) disorder, analyzed in momentum space 
through time-of-flight (TOF) absorption imaging. We ob- 
serve an enhanced density peak in the backscattering di- 
rection (see Fig. [1]), a feature usually attributed to CBS. 
However, we show based on detailed numerical simulations 
that additional effects also lead to an enhanced scattering, 
rendering the direct observation of CBS in our regime a 
subtle issue. In particular, CBS can be a small effect in the 
usual observation channel (momentum space), but appear 
as a well localized bright spot in real space. 

Experiment and simulations. — Our setup [18, en- 
ables the production, in a loffe-Pritchard (IP) magnetic 
trap, of a BEC containing 10^ ^^Rb atoms every 12 s. 
For the experiment described in this paper, a 2D optical 
dipole trap is turned on during the end of the evaporation 
phase. It is made of a vertical standing wave formed by 
crossing two laser beams at a small angle of 4.3° in the 
vertical plane (?/, z). The waists of the employed beams 
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Fig. 1: (Color on line) Transport of a BEC in a 2D disorder. We 
measure the spatial atomic distribution after a 42 ms TOF, for 
various values of the time r spent in the disordered potential. 
Each image, of field of view 1x1 mm, is the average of 30 
shots with different disorder configurations. A peak in the 
backscattering direction is clearly observed for t > 1 ms (as 
emphasized by the circle for r = 1.5 ms). 

define an effective area of 2 x 1.5 mm^ for the 2D waveg- 
uide. The laser wavelength is 766 nm, detuned by 14 nm 
to the blue of the D2 line of ^'^Rb. Once the BEC is 
loaded in the hybrid trap, we decompress the magnetic 
part by lowering the current in the IP trap. This pro- 
duces a rather low density BEC, with a chemical poten- 
tial ji = 270 Hz for N = 10^ atoms. We then release the 
BEC in the 2D waveguide by turning off the IP trap. Due 
to a shift of the magnetic trap center during the decom- 
pression, a velocity along x is communicated to the BEC, 
which can be adjusted by selecting the shut-off time of 
the IP trap. We monitor the subsequent evolution of the 
BEC by absorption imaging. The images can be acquired 
immediately after the 2D trap shut-down to obtain the 
spatial atomic distribution, or after a 42 ms time-of-flight 
(TOF) to measure the velocity distribution. In the con- 
ditions of Figdl the center-of-mass velocity of the initial 
wave packet is 'Uq = 5.9 mm/s (kinetic energy Ek/h = 3.7 
kHz), and the velocity dispersions are Avx = 2.6 mm/s 
and Avy = 1.4 mm/s corresponding to the half- width of 
the Thomas- Fermi distribution. The spatial dimensions of 
the BEC are Ax = 44/im and Ay = 83/im. We study the 
BEC transport in presence of disorder by superimposing 
to the 2D waveguide a speckle pattern, formed by passing 
a vertical beam through a diffusive plate followed by a fo- 
cusing lens. The disorder is here isotropic, contrary to the 
case of [19]. The speckle light is derived from the same 
laser as the 2D trap, frequency shifted by 160 MHz. The 
numerical aperture of the speckle optical system is 0.22, 
resulting in a disorder correlation length ( = 0.54 /im. In 
the plane of the BEC, the speckle intensity envelope is a 
Gaussian of waist = 0.74 mm. The total speckle power is 
89 mW, yielding a speckle potential magnitude Vsp = 2.7 
kHz. The speckle potential is turned on abruptly after 



the release of the BEC, for a variable duration r. Then, 
both speckle and 2D trap beams are switched off and the 
atomic distribution is measured after a 42 ms TOF. 

In order to interpret the experimental results, we have 
performed numerical simulations of the full quantum dy- 
namics of the atomic cloud. The speckle disorder is numer- 
ically generated as explained in [20]. The time-dependent 
Schrodinger (if atom-atom interactions are neglected) or 
the Gross-Pit aevskii (when interactions are important, 
mostly just after the BEC is released from the trap) equa- 
tion is numerically solved using a spatial discretization on 
a grid whose size is much smaller than both the atomic 
de Broglie wavelength and the disorder correlation length. 
Two methods have been used: direct integration of the re- 
sulting coupled time-dependent differential equations, or 
a split technique alternating phases in configuration space 
(where the spatial part of the Hamiltonian is diagonal) 
and momentum space (where the kinetic energy is diag- 
onal). Both methods give the same results. For a given 
realization of the disorder, both the momentum distribu- 
tions and the simulated spatial densities after the TOF 
display strong local fluctuations. This is why we perform 
a configuration averaging over few tens of disorder realiza- 
tions. Since our simulation can either include or exclude 
the effect of interactions in the BEC, we can estimate their 
impact on CBS for the parameters used in the experiment. 
We find that for our parameters interactions are respon- 
sible for a small reduction of at most 10% of the CBS 
contrast. 

Results. — Fig. [T] shows the measured distributions 
for increasing values of r indicated in white. Each image 
actually corresponds to an average of 30 individual shots, 
between each the diffusive plate is rotated to achieve a 
different disorder configuration. The field of view is 1 x 1 
mm. For each image the color scaling is adjusted such 
that it corresponds to the full dynamics of the signal. The 
white arrow on the first image indicates the direction of the 
initial center-of-mass velocity. Qualitatively, we observe 
that the spatial distribution after TOF rapidly develops 
into a ring as r increases, which is consistent with the 
image of an elastic redistribution of the atomic velocities 
due to multiple scattering. As the diffusive component 
(the ring) develops, the initial wave packet centered at 
vq gets depleted. The ring radius is related to vq (we 
measure a radius of 259 /im, consistent with the mean 
atomic velocity) and its width to the velocity dispersion. 
In addition, the abrupt application of the speckle potential 
results in a "disorder broadening" of the ring, of the order 
of 30% of its radius. Within 2 ms the diffusion ring has 
developed fully. However, the most striking observation is 
the clear appearance, after roughly 1 ms, of a peak in the 
backward direction -vq^ the expected signature of the CBS 
effect [H]. 

Fig. [2] shows the numerical results obtained using the 
experimental parameters. The upper panel corresponds to 
the spatial distribution after a 42 ms TOF, and the lower 
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Fig. 2: (Color on line) Numerical results. (A) shows the spa- 
tial distribution after a 42 ms TOF (same field of view as in 
Fig. [T]). (B) is the true velocity distribution (infinite TOF). 
The parameters are those of the experiment (see text). 



one to the true velocity distribution (infinite TOF). The 
spatial and color scales are the same as in Figdl As can be 
seen, the simulation is very similar to the experiment. In 
particular, the dynamics of the population of the diffusion 
circle is well reproduced (see below and Fig. [3lA). Both 
finite and infinite TOF simulations yield an enhanced peak 
in the backscattering direction, but the peak at finite TOF 
is significantly higher and narrower. We will analyze the 
origin of this phenomenon in the Discussion section. 

Fig. [3lA illustrates the dynamics of transport, as mea- 
sured by the ratio of forward to backward scattering. This 
quantity is simply the ratio of the integrated densities of 
the "top half" of the images of Fig.[T]to the "bottom half" , 
as defined by the horizontal line passing through the center 
of the diffusion circle. At times larger than the Boltzmann 
time [2j, one expects this ratio to approach unity. The 
circles correspond to the experimental data, while the line 
is from the numerics using the experimental parameters. 
The shaded area reflects the experimental uncertainty of 
±15% on the speckle intensity. The quantitative agree- 
ment is good. We can extract from the simulations var- 
ious quantities characterizing the transport, such as the 
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Fig. 3: (Color on line) Transport dynamics. We compare in 
(A) the ratio of forward to backward scattering (see text), 
measured experimentally (circles) and numerically (line). The 
shaded area reflects the systematic uncertainty of ±15% on the 
speckle intensity. In B, we plot the angular distribution of the 
radially- integrated density (experimental data of Fig. [J), for 
different values of r (see text for details). The missing data 
around 270° correspond to a saturated area on the images. 
The enhanced backscattering peak is clearly visible for r > 1 



momentum- averaged transport mean-free path i = 7.0 /im 
and transport time rtr = 1.2 ms. We then perform a radial 
integration of the experimental images of Fig. [T] to obtain 
the angular distributions reported in Fig. [3p. For clarity, 
all curves have been scaled to 1 in the forward direction 
(0°). The missing data points around 270° correspond 
to a slightly saturated zone on the images, occuring at a 
fixed position and distorting the data in a limited angular 
range which we thus excluded. The amplitude of the error 
bars shown on the figure corresponds to ±2 standard de- 
viations of the statistical fluctuations. As can be seen, an 
enhanced backscattering peak is significantly observed for 
r > 1 ms. The following section is devoted to the analysis 
of the mechanisms responsible for this enhancement. 

Discussion. — Most of the behaviour observed in 
Fig. [2] is due to the specific properties of the initial wave 
packet created when releasing the BEG from the trap. In- 
deed, as shown in [22], the wave packet roughly keeps the 
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Fig. 4: (Color on line) Contributions to enhanced backscatter- 
ing (numerics). We compare the angular dependency around 
backscattering of the radially-integrated atomic density in four 
different situations: (1) quantum simulation using an initial 
wave packet with the same average velocity vq = 5.9 mm/s 
than in the experiment, isotropic spatial Ax — Ay = 54 fim 
and momentum dispersion Av = 1.3 mm/s and linear position- 
momentum correlation. The propagation time in the speckle is 
1.5 ms; (2) classical simulation with the same initial conditions 
as in (1); (3) same as (2) but without position- momentum cor- 
relation; (4) quantum simulation using the initial state of the 
experiment, as in Fig[2l (A) corresponds to a finite TOP = 
texp, and (B) to the momentum distributions (infinite TOF). 

"inverted parabola" shape of the BEG in the Thomas- 
Fermi regime, but, more importantly, is characterized by a 
strong position-momentum correlation: the fastest atoms 
are located near the edges of the wave packet while the 
slowest ones lie at the center (in the frame co-moving with 
the center of mass of the wave packet). If one neglects the 
small quantum fluctuations (the wave packet is far from 
being minimal), there is even an almost perfect linear re- 
lation between the local velocity v and the distance d from 
the wave packet center. In other words, the initial wave 
packet of spatial size Ax and velocity dispersion Av looks 
as if originating from a point source after a free expan- 
sion time texp = Ax/Av. This linear correlation has two 
important consequences on the backscattered atomic dis- 
tribution after time-of- flight. 

First, it yields a transient enhancement of the atomic 
density in the backscattering direction for TOFs around 
texp- Indeed, if the speckle is applied during a short time 
r such that the atomic displacement inside the speckle 



is negligible compared to Ax, all the particles which are 
backscattered will return to the origin (in the moving 
frame) after a time of flight texp- Despite the initial veloc- 
ity spread, all atoms will reassemble at the same point be- 
cause of the quadratic dispersion relation of matter waves. 
We nickname this effect, reminiscent of the spin echo, the 
backscattering echo (BSE). Contrary to CBS, it is an inco- 
herent effect, maximum when the time of flight duration is 
equal to texp- In the experiment, this time is 17 ms along 
X and 59 ms along y. The 42 ms duration of the TOF lies 
between these values, so BSE is expected to be present as 
well as CBS in our situation. 

Second, the CBS contrast is also expected to be maxi- 
mal for a TOF equal to texp, in analogy with experiments 
in optics where the scattering sample is illuminated with 
a diverging light beam [23 . In this case, the CBS peak 
is observed with maximal contrast at the position of the 
effective point source. The CBS peak observed in the "far 
field" (in our case in momentum space) is the convolu- 
tion of the peak at texp (of angular width :^ {ki)~^) by 
the angular spread of the atomic source Av/vq. To il- 
lustrate these effects, we plot in Fig. H] the results of nu- 
merical simulations corresponding to different conditions: 
(1) a quantum simulation using a correlated initial wave 
packet, with an isotropic velocity dispersion Av = 1.3 
mm/s; (2) a classical simulation using the same correla- 
tions and velocity spread as in (1); (3) same as (2) but 
with uncorrelated position and momentum in the initial 
atomic distribution; (4) the quantum simulation of Figl2] 
using the experimental (anisotropic) momentum distribu- 
tion. We compare the angular distributions of the radially- 
integrated atomic densities after a TOF = tgxp (Fig- ISK) 
and at infinite TOF (momentum distribution. Fig. [3B)- 
Simulations (1) and (4) are expected to include all effects, 
both classical and interferential. Simulation (1) yields for 
TOF = texp a narrow peak due to CBS while the mo- 
mentum distribution exhibits only a broad feature, which 
is consistent with the convolution argument given above. 
This is a new situation where CBS is not visible in the 
usual observation channel (momentum space), but dra- 
matic in real space where it appears as a localized bright 
spot. The characteristic CBS sharp feature of curve (1) 
in FiglHK is broadened in curve (4), due to the "astigma- 
tism" resulting from the anisotropic momentum distribu- 
tion of the initial wave packet. The classical simulation 
with initial correlations (2) yields a triangular peak at 
finite TOF, which is attributed to BSE. Indeed, this fea- 
ture disappears when the correlations of the initial state 
are removed (3). A surprising observation is the resid- 
ual broad peak in the backscattering direction that occurs 
in the classical case without correlations (3). Its width is 
found to be rather independent of the TOF, such that it is 
observed also in the momentum distribution (where BSE 
is absent). We attribute this enhancement to a mecha- 
nism similar to the "opposition effect" [24] , also known as 
"opposition surge" [25] or "shadow hiding effect" - it is 
responsible for the increased brightness of the moon when 
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it is full ^3 - linked to the channelling of atoms inside 
the valleys of the speckle potential. It is a purely classical 
phenomenon which exists only when the average optical 
potential is comparable to the initial kinetic energy, and 
only for relatively short times of the order of few rtr- We 
are currently pursuing a more detailed numerical and the- 
oretical investigation of this unexpected effect. 

Let us summarize the rather complicated situation de- 
scribed above: in the case of a large isotropic veloc- 
ity dispersion Av/vq^ the narrow CBS peak may be ob- 
served by choosing a TOF = texp- BSE, responsible 
for a non-interferential backscattering enhancement, may 
be suppressed by using an initial state without position- 
momentum correlations (such as e.g a thermal cloud). In 
our case of an anisotropic distribution, tgxp is not well 
defined and a sharp CBS peak can not be observed. A 
definite separation from the classical enhancement mech- 
anisms presented above is then delicate. 

Conclusion. — In conclusion, we presented in this pa- 
per our experimental and numerical findings on the quasi- 
2D transport of a BEC in a speckle optical potential, ana- 
lyzed through time-of- flight. We observe a rapid diffusive 
behavior characterized by a transport time of the order 
of 1 ms, accompanied by the apparition of an enhanced 
backscattering peak. These observations are nicely repro- 
duced by our numerical simulations using only experimen- 
tal parameters. These simulations also indicate that the 
observed enhancement of backscattering is both due to 
CBS and two additional classical effects, one of which is a 
transient effect arising from the position-momentum cor- 
relations in the initial atomic wave packet [27]. Further 
improvements in the experiment should allow the quanti- 
tative study of the dynamics of CBS and the determina- 
tion of the coherence length, a key parameter in localiza- 
tion experiments. The present work deals with a rather 
dilute BEC, but the chemical potential can be increased 
by more than one order of magnitude to study the in- 
terplay of weak localization and non linearities. Further 
studies could include the search for the "coherent forward 
scattering" effect predicted at the onset of the Anderson 
localization regime [28] . 
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